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ABSTRACT 

We introduce an alternative hypothesis to explain the very low luminosity of 
the cool (L-type) companion to the ~25 Mj up ~8 Myr-old brown dwarf 2M1207A. 
Recently, Mohanty et al. (2007) found that effective temperature estimates for 
2M1207B (1600 ±100 K) are grossly inconsistent with its lying on the same 
isochrone as the primary, being a factor of ~10 underluminous at all bands 
between I (0.8 /im) and V (3.6 /im). Mohanty et al. explain this discrepency 
by suggesting that 2M1207B is an 8 Mj up object surrounded by an edge-on disk 
comprised of large dust grains producing 2.5 m of achromatic extinction. We offer 
an alternative explanation: the apparent flux reflects the actual source luminosity. 
Given the temperature, we infer a small radius (~49,000 km), and for a range of 
plausible densities, we estimate a mass < Mj up . We suggest that 2M1207B is a 
hot protoplanet collision afterglow and show that the radiative timescale for such 
an object is > ~1% the age of the system. If our hypothesis is correct, the surface 
gravity of 2M1207B should be an order of magnitude lower than predicted by 
Mohanty et al. (2007). 



Subject headings: circumstellar matter — planetary systems : formation — plan- 
etary systems : protoplanetary disks — stars: individual (2MASSW J1207334- 
393254) — stars: low-mass, brown dwarfs — stars: pre-main-sequence 
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Introduction 



While radial velocity surveys and other techniques have yielded over 200 extrasolar 
planets, there are to date no convincing images of an extrasolar planet around a star. A few 
candidate "planetary mass objects" have been identified co-moving with p re-main sequence 
stars and young brown dwarfs (e.g. GQ Lup, Oph 162225-240515, etc.; e.g. iNeuhauser et al. 



20051 ; iJayawardhana fc Iyanovll2006r). However further observations have shown these objects 
to be high er mass (e.g. iLuhman et al.l l2007al ). A possible exception is the companion to 



2M1207A f Chauvin et al. 



2004 



2MASSW J1207334-393254 = TWA 26), for which mass 
estimates have ranged between ~3-8 Mj up . Although its mass inferred from evolutionary 
models is well below the deuterium-burning limit, it is unlikely that the object could have 
formed through classical core accretion in a circumstellar disk. Given its estimated mass 
ratio (~ 0.2-0.3), and separation (~ 50 AU) it is thought to have formed via gravitational 
fragmentation, similar to a binary star system (ILodato et al.ll2005l ). 



The properties of the 2M1207 system are discussed in detail in iMohanty et al.l (120071 ). 
The most striking result from that study is that while both the colors and near-IR spectrum 
of 2M1207B are consistent with models of an unreddened ~1600K object, the inferred 
luminosity for the object is ~2.5 mag (factor of ~10) below that expected for a ~5-10 
Myr object at all wavelengths. Even more remarkable, the object is slightly fainter than the 
observed sequence of K a n d L' absolute m a gnitu des for older field objects with T eff ~ 1600 K 
( IGolimowski et al.l 120041 ) . IMohanty et al.l (120071 ) rule out a handful of simple resolutions to 
explain the apparent underluminosity of 2M1207B, and ultimately settle for an unlikely but 
testable hypothesis: that the object is obscured by an edge-on disk of large circumstellar dust 
grains producing 2.5 m of gray extinction. In this contribution, we propose an alternative 
explanation for the low luminosity of 2M1207B: namely that it has a small radius. In 
this scenario 2 M1207B is the ho t, long-lived afterglow of a recent collision between two 
protoplanets (cf. IStevensonlll987i). The observational consequences of pr o toplanet collision s 
have been discussed elsewhere ( jSternl Il994l ; IZhang fc Sigurdssonl 120031 ; lAnic et al.l 120071 ). 
Here we review the problem of 2M1207B, introduce a new hypothesis, discuss its merits and 
deficiencies, and offer an observational test that can rule it out. 



2. The Luminosity of 2M1207B 



We first review some observational properties of the 2M1207 system. 2M1207A i s 
classified as M8 with emission-line activity characteristic of T Tauri stars (IGizisI 120021 ). 
The syste m appears to harbor a circumstellar accretion disk as e videnced by broad Ho; 
emission (IMohanty et al.l 120031 ) . mid-IR excess (jSterzik et al.l 12004 ). and outflow activity 
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( IWhelan et al.l 120071 ). The motion of 2M1207A is consistent with membership in the TW 
Hya Association (TWA), a l oose group of ~20 stars with mea n age ~8-Myr-old situated at 
a mean distance of ~50pc (IWebb et al.lll999l ; iMamajeki 120051 ) . The optical/near-IR colors 
of 2M1207A are consistent with no reddening given its spectral type. Comparison of its 
position in the H-R diagram with theoretica l models suggests a m ass of ~25 Mj up and age 
consistent with its membership in the TWA (IMohanty et al.l 120071 ) . 



2M1207B shares common proper motion and parallax with 2M1207A within the as 



trometric uncertainties (mil l iarcsecond- level mea surements over a few years; IChauvin et al. 



2004 120051 : ISong et al.ll2006l : Mohantv et al.ll2007[ ) . The idea that B could be a foreground or 
background field L dwarf has been ruled out. At a dista nced of 66 ± 5 pc , the c ompanion is 
at a projected separation of 51 ±4 AU from the primary. IMohanty et al.l (120071 ) concentrate 
their analysis of the temperature for 2M1207B on comparison of data to the Lyon group 
model atmospheres. From the available H and K-band spectra, they d etermine the best fi t 
model atmosphere has T cfr = 1600 ± 100 K from the DUSTY grid of kllard etaD (boOlh . 
A range of surface gravity from log(g) = [3.5-4.5] was explored, but the value is poorly 
constrained. They also compare the available photometry of 2M1207B from 0.9-4.0 /im 
with predictions from the DU STY models, a nd co nclude that they are consistent with this 
temperature estimate. While iLeggett et al.l (120011 ) demonstrate th a t DU STY models are a 
somewhat poor fit for old late L-type field dwarfs, IMohanty et al.l (120071 ) have shown that 
they produce adequate spectral fits for the young objects 2M1207B and AB Pic B. 2M1207B 
is significantly redder and dustier than typical L dwarfs, as predicted for low surface gravity 
objects. 



We explore a complementary approach, comparing the spectra published in IMohanty et al 



(120071 ) as well as the available photometry with template objects drawn from wide field 
surveys (predominantly older objects). The low resolution H and K-band spectrum avail- 
able at signal-to-noise ratio of 3-10 is morphologically similar to other L dwarfs suspecte d 



of having low gravity: 2 M ASS J01415823 - 46335 74 (2M J0141; iKirkpatrick et all 12009 1. 
2MAS S J22244381-015852 1 JCushing et al.l l2005h . and SDSS J22443167+2043433 (SDSS 



Allers et al. 



2007 



(jCushing et al. 



J2244; iKnapp et al.l I2004T). AU of th ese objects exhibit weak metal resonance lines (e.g. 



2005 



McLean et al 



Gorlova et al.l 12 003). stronger than expected CO for their spectral type 



20031 ). and unusual pseudo-continua in the H-band spec- 



1 Our adopte d distance to 2M1207 via the cluster parallax method has increased b y 2a compared t o 



Mamaiekl (|2005l . 53 ±6 pc) due to the effects of a improved proper motion for 2M1207 ijSong et all 12006). 
and a revised estimate of the TWA group velocity (22.4 ± 1 kms -1 ; which follows Mamajek 2005, but omits 
th e deviant paralla x for TWA 9). The updated velocity increases the distances to the other TWA members 



Mamaiekl (|2005l) by 7%. 
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tra attributed to collision induced molecular hydrogen absorption (e.g. iBorysow et al.lll997l ). 
all indications of low surface gravity. Given the morphological correspondence between the 
spectrophotometry of 2M1207B and these other low gravity L dwarfs, it is reasonable to 
assume that 2M1207B has a similar nature. Further, the spectrum of 2M1207B shows no 
signs of CH4 absorption which would indicate T e ff below 1400 K. 



The colors of 2M12 7B are very red compar ed to observed sequences of field L dwarfs 



le.g. 



Knapp et al.l 120041 ; iGolimowski et al.l 120041 ) . The anomalous L dwarfs listed above 



also exhibit this b ehavior which can be attributed to low gravity (e.g. iBurrows et al.l 12006 



Allard et al.ll200ll ). Taking the observed (J-H) colors of 2M1207B and several plausible 
intrinsic colors matches, we searched for re ddening solutions that would fit the SED of 
2M1207B. Adopting the colors of 2M J0141 (IKirkpatrick et al-lhoOO ) as a low gravity early 
L template, we derived Ay ~ 9 mag, which matches t he JHK photometry well. However, 
2M J0141 has T cfr = 2000 K fcirkpatrick et alil2006h . and A v ~ 9 mag implies A K ~ 1 
mag, which would be insufficient to move 2M1207B ab ove the old dwarf s equence for that 
T e ff. Alternatively, adopting the colors of SDSS J2244 (IKnapp et al.l 120041 ) as a late-L low- 
gravity template, we arrive at Ay ~ 3.7, again reproducing the colors of 2M1207B within 
the (rather large) errors from ~l-4 /xm. However, the required reddening (A# ~ 0.4 mag) 
is insufficient to solve the underluminosity of 2M1207B. If we assume that 2M1207B has the 
{J-H) colors of a late-type M dwarf, we derive Ay — 11 mag, but cannot reconcile the 
observed colors without invoking excess emission in the K and L-bands. In t he limit of zero 



extin ction, one can find models of extremely low gravity that fit the SED (e.g. lMohanty et al. 
20071 ). so we take that as the simplest assumption consistent with the observed properties of 
known L dwarfs and informed by model atmospheres. 

In order to estimate the bolometric luminosity of 2M1207B, we must also estimate an 
appropriate bolometric correction to apply to the observed absolute magnitude. Given the 
distance to the source, the lack of evidence for interstellar reddening, the available photome- 
try, and the temperature estimate discussed above, we can a pply a bolometric correction to 
any flux estimate from 0.9-4.0 //m. IGolimowski et al.l (120041 ) demonstrates that BCj^ varies 
little as a function of spectral type for L dwarfs, so we apply the X-band BC to 2M1207B 
to minimize th e uncertainties in the estimate of (BCk = 3.25 ±0.14 mag; same as in 
Mamajekll2005l ). This results in a luminosity estimate of log(L/L ) = -4.54 ± 0.10 dex. For 
T eff = 1600 ± 100 K, the DUSTY models predict BC^ = 3.56 ± 0.07 mag, and a luminosity 
of log(L/L Q ) = -4.66 ± 0.08 dex. The difference between adopting the empirical or the- 
oretical BC values is within the errors, so we conservatively adopt the empirically derived 
log(L/L ) as an upper limit. Given its luminosity, 2M1207B lies 4-7 times below where 
ex pected for an age range of 5-1 Myr given its inferred temperature range of 1500-1700 
K. iMetchev fc Hillenbrand! (120061 ) derive an effective temperature for the young (~0.3 Gyr) 
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L/T transition object HD 203030B that is considerably lower (by ~230 K) than predicted 
for its spectral type (L7.5). If 2M1207B were in fact a late-L spectral type with a tempera- 
ture of 1100 K, we could reconcile its position in the H-R diagram given its age. However, 
this is ^500 K cooler t han the T e ff derive d from spectral synthesis models according to 



Mohanty et al.l (120071 ). iBurrows et al.l (120061 ) suggest that the temperature of L-T transition 
objects should only weakly depend on temperature. Exploring this solution to the 2M1207B 
problem requires higher SNR spectra and further study. 

Can the theoretic al evolutionary t racks be wrong in luminosity by such a large fac- 
tor? As pointed out by lMohanty et al.l (120071 ). another well-studied young, low- mass binary 
(AB Pic) exh i bits H RD positions consistent with the age of the group (Tuc-Hor; ~30 Myr). 
Marley et al.l (120071 ) have argued that the luminosities of young planets formed through core- 
accretion are likely to be over-predicted in models that initially start objects in high entropy 
states. However, this argument does not apply if the obj ect for med through gravita tional 
fragmentation, as has been suggested (ILodato et al.ll2005l ). The iMarley et al.l (120071 ) mod- 
els predict that for ages >Myr after accretion has ceased, all planets that formed through 
core-accretion with masses of <10 Mj up will be colder than T e fj ~ 800 K and with log(L/L ) 
< -5.2 dex. Hence, while 2M120 7B is vastly underlum i nous for its T e ff com pared to the 
"hot-start" evolutionary models of IChabrier et al.l (120001 ); iBaraffe et al.l (120031 ). it is overlu- 
minous and too hot for the fiducial "cold-start" core accretion models for <10Mj up from 
Marley et al.l (120071 ). Observationally, there seems to be no trend that would suggest an 
error in the evolutionary tracks that could account for the HRD position of 2M1207B. 



As pointed out in iMohanty et al.l (120071 ) one cannot reconcile the observed spectrum 
and SED of 2M1207B with its apparent low luminosity given the available models. Here we 
take a different approach, adopting the derived temperature of 1600 K, and postulate that 
the observed source flux is a reflection of its actual luminosity. 



3. A Planet-Collision Theory for 2M1207B 

We hy pothesize that 2M1207B is the result of a recent collision of two protoplanets (cf. 
Sterol 1 19941 ). Adopting the temperature and luminosity in §2, one derives a radius of 48,700 
± 8,800 km (= 0.68 ±0.12 Rj up , 7.6 ±1.4 R e ). For a range of densities, the inferred mass 
and gravity are given in Table [TJ 



We hypothesize that the object is a hot protoplanet of density ~1 gem 3 . With this 
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radius and densitj]^, the fiducial planet would have a mass of Mb = 81 M ffi (0.85 Ms ainrn = 4.7 
M Neptune ). Throughout this discussion we will refer to hypothetical protoplanets B 1 and B 2 
which merged to produce body B (2M1207B). In our solar system, the largest planetesimals 
to impact the planets during the late stages of accretion appear to have had mass ratios of 
order 7 = MbiIMbx ~ 0.1, in cluding the bodies respons ible for producing the obliquities of 
Saturn, Uranus, and Neptune (ILissauer fc Safronovlll99ll ). Thus we assume that our ~81 M® 
object was the product of a collision between protoplanets with masses Mm = (1 + 7)~ 1 Mr 
and M B2 = 7(1 + 7)~ 1 M B , or M B1 = 74 M e and M B2 = 7 M ffi . Followinglwetherilj fjl98oh . 



the minimum impact velocity for two planets will be their mutual escape velocity, defined 
as: 



2G(M B1 + M B2 ) 



J mut 



Rbx + Rb2 

In the simplified case of identical densities p for bodies Bi, B 2 , and B, one can simplify this 
equation in terms of mass and radius of the final planet B: 



Vmut 



11.2 km s" 



M, 



B 



M a 



Rb 

Rffi 



(1 + 7)* 
, 1,1 
(1 + 73)2 



(2) 



where the last factor is within <15% of unity for all 7 < 10 . In our fiducial model, the radii 
of the fiducial impactors are 3.0 and 6.4 R m , respectively. This leads to a fiducial impact 



velocity > 30.6 kms l . Following ISternl (119941 ). one can calculate the radiative timescale of 
a long-lived afterglow from the collision of bodies B\ and B 2 . 



Trad ~ 0.33 Myr 



M 



B2 



Rb 



1000 u io 



(3) 



where V\q is the impact velocity of the impactor in units of 10 kms -1 . Rb is the radius of 
the planet after collision, and is the radius of Earth. T 100 o is the temperature of the 
emitting photosphere in units of 1000 K. We can rewrite the radiative timescale in terms of 
the impactor mass ratio 7, the properties of the final body B, and assuming impact velocity 
equals v mut : 



r rad ~ 0.41 Myr 



Mi 

Ma 



Rb 

R(f, 



rp — 4 



1000 



/ 



(4) 



where / = 7(1+7) 3(1 + 73) 1 . For 7 <C 1, / ~ 7 (to within <40% accuracy for 7 < 10 1 ). 
For our fiducial model, r rad > 59 kyr or ~1% the age of the TW Hya association. These 



2 The giant planets in our solar system have bulk densities of ~ 0.7-1.4 gem 3 , and the known transiting 
hot Jupiters have d ensities of ~0. 3 -1.3 g cm~ 3 (jBakos et alJl20071 ). The post-accretion evolutionary track of 
a 1 M Jup object bv lMarlev et al.l (|2007h has density -0.5-0.6 gem" 3 in its first 10 Myr. 
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radiative timescales are not negligible, and suggest that a hot afterglow could be visible for 
an appreciable fraction of the system lifetime. I n mo deling the collision of a Jupiter and 
an Earth-like protoplanet, IZhang fc Sigurdssonl (120031 ) estimate that less than 1% of the 
impact energy is radiated away in the initial prompt flash. They further argue that most 
of the collision-deposited energy is locked up deep in the post-collision planet, and radiated 
over a long timescale as an afterglow that peaks in the IR. We propose that 2M1207B is 
such a long-lived afterglow. 

Could a plausible circum(sub) stellar disk form a planetary system with the required 
properties? The formation and collision of two such large planetesimals at radii > 10 AU 
in a protoplanetary disk surrounding a b rown dwarf is very un likely given the mass surface 
density and orbital timescales expected (jGoldreich et al.l 12004 ). Perhaps 2M1207B formed 
at smaller radii as the ice-line i n the disk of 2M1207A swept throu gh a large range of 
inner radii from 10 to 0.1 AU (cf. Kennedy. Kenyon. fc Bromley! 120061 ) as the young brown 
dwarf evolved (see iBossI 120061 . for an alternate scenario). If we consider a primordial disk 
surrounding 2M1207A that is marginally gravitationally stable ( M ^ k ~ 0.1) it would have 
a total ga s +dus t mass of 2-3 Mj up . Adopting the protoplanetary core mass scenario of 
Ida fc Linl ([2004], see also Lodato et al. 2005), the time evolution of the mass of a planet 
accreting 10 18 g planetesimals is: 

3 / v \ 21 / 5 

M p (t) » 8M S ' 



10 6 yr 



a 



10 g cm 2 

-9/5 f M A 



1/2 



(5) 



1AU/ \M Q/ 

Where t is time, is disk surface density of solids, a is the orbital distance, and Ma is 
the mass of 2M1207A. Assuming a disk with the above mass, mass surface density profile 
£ oc a" 1 , and an outer radius of 15 AU, the total disk mass surface density at 3 AU is 
250 gcm~ 2 . If we consider a gas to dust+plus ice ratio of 25 (100/4), we arrive at a mass 
surface density in solids of 10 gem -2 at 3 AU. Using the above equation for a brown dwarf 
of mass 0.025 M©, we estimate that a core of 5-10 M e can form within ~3 Myr at this 
radius. In 10 Myr, a similar mass core could form at a distance of 5 AU from the brown 
dwarf. Assuming the above disk model parameters, roughly half of the total disk mass (solid 
and gas) resides inside of 7.5 AU and half outside. It is at least plausible that 2-4 cores of 
5-10 M e could form between 1-10 AU within 3-10 Myr in this system utilizing the bulk of 
available solids in the system. If two of those cores accreted enough gas to form Neptune- 
to-Saturn mass protoplanets, we can envision a scenario where: a) failed cores of 5-10 M e 
could collide with a successfully formed gas/ice giant protoplanet, creating the observed hot 
collisional afterglow; and b) another gas/ice giant, along with the presence of the remnant 
primordial disk, could eject 2M1207B to its observed orbital radius of 50 AU. Motivated by 
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evidence for gas giants at large separations having created observed structure in debris disks, 



Veras fc Armitagd (120041) have investigated gas giant migration/ejection scenarios in disks. 



Thommes et al.l (120031 ) have also proposed that Neptune and Uranus formed closer to the 



Sun (between Jupiter and Saturn) and were ejected to larger orbital radii through dynamical 
processes. We note that the remnant disk s urrounding 2M1207A has a mass comparable to 
that we propose for the ejected 2M1207B (IRiaz fc Gizisl 120071 ) though its outer radius is 
unconstrained from current observations. This is, of course, a highly improbable series of 
events. 



4. Predictions 



The hypotheses of whether or not 2M1207B is a hot protoplanet collision afterglow or 
is obscured by a dense disk of large dust grains can be tested. In the scenario proposed here, 
2M1207B is actually a ~80M e object with radius ~49,000km. The surface gravity of such 
an object in cgs units would be log(g) ~ 3 (Table 1). Thi s is significantly lower than than for 
5-10 Myr-old 3-8 M Jup objects which have \og(g) w 4 (IMohantv et al.lboOTl ). If 2M1207B 
possesses an edge-on disk exhibiting grey extinction, then spatially resolved ground-based 
observations should reveal: a) an infrared excess at A > 4 fim from the disk; b) a 10 jttm 
silicate absorption feature consistent with the disk being edge-on; c) polarized emission from 
scattered light at shorter wavelengths; and/or d) resolved sc attered light emission consistent 
with an edge-on dust disk system (cf. iLuhman et al.ll2007bl ). 



If 2M1207B is actually a physically smaller (and therefore lower mass) companion, it 
sh ould exhibit nea r -infra red spectra: a) consistent with the 1600 K temperature advocated 
by IMohantv et al.l (120071 ): and b) low surface gravity (log(g) ~ 3) in high S/N spectra. As 
mentioned above, surface gravity affects the spectra of very c ool objects in ways th at can be 



observed through analysis of atomic and molecular features. iGorlova et al.l (120031 . see also 



Kirkpatrick et al. 2006) suggest that log(g) ca n be estimat e d to w ithin 0.3-0.5 dex from high 
S/N near- infrared spectra of M and L dwarfs. lAllers et al.l (120071 ) specifically investigate the 
gravity dependence of the Na I feature at 1.14 /iin while Gorlova et a l. provide a preliminar y 
calibration of the surface gravity effects of K I at 1.25 ixm (see also iMcGovern et al.l I2004T ) . 
These effects should be clear in modest S/ N spectra (20-30) easily distinguishing between the 
log(g) ~ 4 model of lMohanty et al.l (120071 ) and the log(g) ~ 3 model proposed here. Further, 
we anticipate that our protoplanetary collisi on remnan t would be metal-rich compared to the 
primary. Models from lBurrows et al.l (120061 ) as well as lFortney et al.l (120061 ) demonstrate the 
significant differences in brown dwarf and gas giant planet atmospheric models by varying the 
metallicity. Such effects would be easily observable in S/N ~ 20-30 spectra, obtainable in 1-2 
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nights of observing time on a 6-10 meter telescope equipped with adaptive optics. Perhaps 
future surveys will uncover additional hot protoplanet collision afterglow candidates with 
even smaller inferred masses. 
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Table 1. Predicted Quantities for 2M1207B 



p Mass Mass log(<?) 
(gem' 3 ) (M ffi ) (M Jup ) (cms" 2 ) 



0.5 


40.5 


0.13 


2.83 


1.0 


81.0 


0.25 


3.13 


1.5 


121.5 


0.38 


3.31 


2.0 


162.0 


0.51 


3.43 



